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Abstract 

We study the long time behavior of an Ornstein-Uhlenbeck process under the influence 
of a periodic drift. We prove that, under the standard diffusive rescaling, the law of the 
particle position converges weakly to the law of a Brownian motion whose covariance 
can be expressed in terms of the solution of a Poisson equation. We also derive upper 
bounds on the convergence rate. 

1 Introduction 

In this paper we study the long time behavior of solutions of the following Langevin 
equation: 

T£{t) = v{x(t))-x(t) + a${t) , x{t)€W, (1.1) 

where f3(t) is a standard Brownian motion and a, t > 0. The parameter r can be 
thought of as a nondimensional particle relaxation time, which measures the inertia 
of the particle. The drift term v is taken to be smooth, periodic with period 1 in all 
directions; further, it is assumed that it satisfies an appropriate centering condition. 

It is well known that as r tends to the solution of converges with probability 
1 to the solution of the Smoluchowski equation 

z(t) = v(z(t)) + a/3(t) , x(t) e R" , (1.2) 

uniformly over every finite time interval, see e.g. INel67l Ch 10]. The problem of ho- 
mogenization for equation ( ll.2l i has been studied extensively over the last three decades 
for periodic I BLP78 Bat85 Par99 1 as well as random | CX97 KOOl LOY98| drifts. 
For the case where v(z) is a smooth, periodic field which is centered with respect to the 
invariant measure of the process, it is not hard to prove |BLP78 Ch 3] that the rescaled 
process ezit/e^) converges, as e tends to 0, to a Brownian motion with a positive def- 
inite covariance matrix K,. The proof of this functional central limit theorem is based 
on estimates on the spectral gap of the generator of the process z{t). 

The long time behavior of particles with non-negligible inertia, whose evolution is 
governed by equation M.W has been investigated by Freidlin and coworkers in a series 
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of papers llhW!j8l IPreOll IFWOTI IFW 991. Among other things, Hamiltonian systems 
under weak deterministic and random perturbations were studied in these papers: 

Tx = ~\/V(x) + e(-Kx + -f) + y/ecr(3, (1.3) 

with K, 7 e M. It was shown that, under appropriate assumptions on the potential 
V(x), the rescaled process {x{t/e), y(t/e)} converges weakly to a diffusion process 
on a graph corresponding to the Hamiltonian of the system H = ^tx^ + V{x). 

On the other hand, the problem of homogenization for M.W has been investigated 
less. This is not surprising since the hypoellipticity of the generator of the process M.W 
renders the derivation of the necessary spectral gap estimates more difficult. Homog- 
enization results for the solution x{t) of M.W have been obtained, to our knowledge, 
only for the case where the drift v{x) is the gradient of a potential. In this case the in- 
variant measure of the process {x(t)^ x{t)} is explicitly known and this fact simplifies 
considerably the analysis. This problem was analyzed for periodic IRod89l as well as 
random potentials fPV85l. In both cases it was shown that the particle position con- 
verges, under the diffusive rescaling, to a Brownian motion with a positive covariance 
matrix JC. The proofs of these homogenization theorems are based on the techniques 
developed for the study of central limit theorems for additive functionals of Markov 
processes IKV86I, together with a regularization procedure for appopriate degenerate 
Poisson equations. Related questions for subelliptic diffusions have also been investi- 
gated rNor94I INor97llBBJR95l . 

The purpose of this paper is to prove a central limit theorem for the solution of 
the Langevin equation M.W with a general periodic smooth drift v{x) and, further, to 
obtain bounds on the convergence rate. The proof of our homogenization theorem relies 
on the strong ergodic properties of hypoelliptic diffusions. The techniques developed 
in IEPRB99 EHOO | enable us to prove the existence of a unique, smooth invariant 
measure for M.W and to obtain precise estimates on the solution of the Poisson equation 
— £/ — g, where C is the generator of the process M.W and the function g is smooth 
and centered with respect to the invariant measure. Based on these estimates it is rather 
straightforward to show that the rescaled particle position exit/e^) convergences to a 
Brownian motion, using the techniques developed in |KV86|. Obtaining bounds on 
the rate of convergence requires more work. For this, we need to identify the limiting 
Brownian motion and to introduce an additional Poisson equation. 

The sequel of this paper is organized as follows. In section |2] we introduce the 
notation that we will be using throughout the paper and we present our main result, 
Theorem l2.1l In section|3lwe prove various estimates on the invariant measure of M.W 
and the solution of the cell problem, and we also derive estimates on moments of the 
particle velocity. The proof of Theorem l2.1l is presented in section]?] Finally, section|5l 
is reserved for a few concluding remarks. 
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2 Notation and Results 

Consider the following Langevin equation in M" : 

Tit{t) = v{x{t)) - x(t) + al3(t) , (2.1) 

with initial conditions a;(0) = x, x{0) = (y^)^^j/. We assume throughout this pa- 
per that V S C°°(T"). Introducing y{t) — ^/rxit), we rewrite ( 12. 1> as a first order 
stochastic differential equation: 



dx{t) — —=y{t) dt , 

dy(t) = -^v{x(t)) dt - -y(t) dt + -^d(3{t) 
\Jt t 

We denote by £ the generator of the process {x{t), y(t)}: 



(2.2) 



1 , _ , , _ . 1 / _ cr 



^2 



£ = (y • V, + v(x) • Vy) + - \^~y ■ Vy + y^v) ' ^^'^^ 

By Theorem l3 . 1 I below. the process {x{t), yit)} admits a unique, smooth invariant mea- 
sure, denoted by ii{dx, dy). 

Consider now the cell problem 

-£$ - . (2.4) 

This equation has a unique, smooth solution in the appropriate function space by Theo- 
rem]^], provided that J v{x) fiidx^ dy) — 0. We define the symmetric, positive n x n 
matrix /C such that 



IC-' = — Vy^®Vy<i>d^i. (2.5) 



r 

The main result of this paper is that the particle position, under the standard diffusive 
rescaling, converges weakly to a Brownian motion with covariance JC^. We further- 
more give upper bounds on the rate of convergence in the following metric. Let B 
denote a separable Banach space and B* be its dual space. Given two measures /ii and 
/i2 on B, we also denote by C{pi, ^2) the set of all measures on B^ with marginals /ii 
and /i2. With these notations, we define the following metric on the space of probability 
measures on B with finite p-moment: 

P . r f \i(x)~£(y)\P ^ , , ^ 

lAii -M2ip = sup mf / — — ne(dx,dy). (2.6) 

This distance is close in spirit to the p-Wasserstein distance 



IIImi -M2|||pM/ = J \\x - y\\P fiidx,dy) , 

so we will refer to it as the weak p-Wasserstein distance. Note that the distance \2.6\ 
gives a locally uniform bound on the distance between characteristic functions x^(^) — 
J e*^<^' ^i(dx): 

IXmi(^)-Xm.(^)I< IKII lllMi-M2|||p . 
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In particular one has |||/ii — /i2 = if and only if /ii = /i2- 

In order to simplify notations, we define the fast processes = y(e~^t) and 
= x(e^^t). We will also from now on use the notation B = C([0,r],R"), for 

a value T > that remains fixed throughout this paper. Now we are ready to state the 

homogenization theorem. 

Theorem 2.1 Let x{t) be the solution of M.W . in which the velocity field v £ C°°(T") 
satisfies J v{x) fx(dx, dy) = 0. For T > fixed, denote by ii^ the measure on B given 
by the law of the rescaled process ex^ and by fi the law of a Brownian motion on R" 
with covariance K? as defined in \2.5\ . Then, for every p> 1 and a G (O, 5), there is 
a constant C such that 

- Mlllp < Ce^ , (2.7) 

for all e e (0, 1). Furthermore, if one denotes foy tt^ : B ^ C([0, T], R) the projection 
given by {'Kkx){t) = {k, x{t)), one has the bound 

IlkfeMe - T^*M\p,w ^ C'e" . (2.8) 
for every fc G M" with \\k\\ < 1. 

Remark 2.2 The condition / v{x) ^{dx, dy) = ensures that there is no ballistic mo- 
tion involved. In the general case, one can write v — J v{x) ii{dx, dy) and define 
eXf — ex{e^^t) — e^^vt. Then, Theorem l2. ll holds for exf. 

Remark 2.3 If n = 1, the bound ( 12. 8t is much stronger than the bound il.li . If ?i > 1 
however, this bound does not imply any form of convergence /i^ ^ /i. It is indeed 
possible to construct two Gaussian stochastic processes x{t) and y(t) with values in R^ 
such that the law of x differs from the law of y and such that, for every k E R'^, the law 
of (fc, x) is identical to the law of (fc, y). As an example, choose three i.i.d. Gaussian 
centered random variables ai , 02 , 03 and define 

Xiiti) = ai X2{ti) = a2 Xi(t2) = 03 2:2(^2) = Oi 
yi(ii) = fli 2/2(^1) = 0.2 yiih) = 02 y2(t2) = 03 ■ 

It is an easy exercise to check that these two processes possess the required properties. 

Remark 2.4 Convergence in the weak p-Wasserstein distance alone does not imply 
weak convergence, as the space of probability measures on B is not complete under 
III • lip. This can be seen by taking B = £^ and choosing for /i„ the Gaussian measure 
with covariance 

g„==diag(l,i,...,i 0,...) . 

It is straightforward to check that this forms a Cauchy sequence with respect to ||| • 1^, 
but does not converge to any measure supported in i'^. (It does however converge 
weakly to a limiting measure in a weaker topology, and this is always the case.) In our 
case, it is easy to check that the sequence of measures /i^ is tight, since the generalized 
Kolmogorov criteria IRY99I Thm 2.1] provides us with uniform bounds on the a- 
Holder constant (with a < i) of the process x^. Tightness, together with convergence 
in the weak p-Wasserstein distance then implies weak convergence. Note also that 
even though convergence in the weak p-Wasserstein distance alone does not imply 
weak convergence, it does imply weak convergence of finite-dimensional distributions. 
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Remark 2.5 The covariance, or effective dffusivity, K? of the Hmiting Brownian mo- 
tion depends on the a and r. It is shown in fPS03l that as t tends to the covariance 
converges to the one obtained from the homogenization of equation 1 1.21 We re- 
fer to I PS03 1 for further properties of the effective diffusivity, together with numerical 
experiments for various fields v{x). 

Remark 2.6 For simplicity, we choose the molecular diffusion cr to be a constant 
scalar Taking for a a positive definite matrix would only require a slight change in 
our notations. We could even allow a to depend on a; in a smooth way, as long as it 
remams strictly positive definite for all x e T". The results from |EPRB99 EHOO| 
then still apply and one can check that all the bounds obtained in section Instill hold. 
Since the proof of Theorem 12. II itself never uses the fact that a is constant, all of our 
result immediately carry over to this case. 

Remark 2.7 For simplicity, we assumed the initial condition {x^y) to be determin- 
istic. However, it is easy to check that all our arguments work for randomly dis- 
tributed initial conditions provided that they are independent of the driving noise and 
that E exp (5|| < oo for all 5 G (0,a~^). In particular, one can take the initial 
condition to be distributed according to the invariant measure. 

The proof of this theorem will be presented in section|4] 

3 Preliminary Estimates 

In this section we collect various estimates which are necessary for the proof of the 
homogenization theorem. In section lTTI we study the structure of the invariant measure 
fi for (12. 1> . We show that it possesses a smooth density with respect to the Lebesgue 
measure and we derive sharp bounds for it. Further, we investigate the solvability of 
the Poisson equation 

-Cf = h, (3.1) 

where his a smooth function of x and y which is centered with respect to /u. We prove 
that equation (13. 1> has a smooth solution which is unique in the class of functions 
which do not grow too fast at infinity. 

In section l3^ we derive estimates on exponential moments of the particle velocity. 
Roughly speaking, these estimates imply that the particle velocity grows very slowly 
with time. 

3.1 Bounds on the invariant measure and on the solution of the Poisson equation 

If w = 0, the invariant measure for i2.l\ is given by fi = "'^ dx dy. This is "almost" 
true also in the case w ^ 0, as can be seen by the following result. 

Theorem 3.1 Let /i be the invariant measure for H2.ll and denote by p(x, y) its density 
with respect to the Lebesgue measure. Then, for every 5 € (0, 2a~^) one can write 

p(x,y) = e-^\\y\\"g(x,y) , geS, (3.2) 
where S denotes the Schwartz space of smooth functions with fast decay. 
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Proof. The proof follows the lines of [EPRB99 EHOOl. Denote by 0t the (random) 
flow generated by the solutions to M.W and by Vt the semigroup defined on finite 
measures by 

(7'tAi)(^) = E(/io0,-i)(^). 

Since 9t + /! is hypoelliptic, Vt maps every measure into a measure with a smooth den- 
sity with respect to the Lebesgue measure. It can therefore be restricted to a positivity 
preserving contraction semigroup on L^(T" x M", dx dy). The generator L of Vt is 
given by the formal adjoint of C defined in ( 12. 3> . 

We now define an operator K on L^(T" x M", dx dy) by closing the operator de- 
fined on by 

X = £e-|llyf . (3.3) 

The operator K is then given by 

5(J^\.. „o 1 



K = -|-A, + ^(1 - ^) hf + HScr' -l){y V, 



1 71 

+ ^(yV, +f(2;)-V.y)-— . 

Note at this point that S < 2a^^ is required to make the coefficient of ||y||2 in this 
expression strictly positive. This can be written in Hormander's "sum of squares" form 
as 



Xi = dy if i — 1 ... n, 



with 



^« = V T ( 1 - — ) U^-n if i = (n + 1) . . . 2n, 

^0 = -{Sa' -l)[y V.y + - j + -^(2/ • V, + v(x) -Vy)- — . 

Since t; is C°° on the torus, it can be checked in a very straightforward way that the 
assumptions of |EHOO" Thm. 5.5] are satisfied with = I - - Ay + \\y\\^. 
Combining this with lEHOOl Lem. 5.6], we see that there exists a > such that, for 
every 7 > 0, there exista a positive constant C such that 

IIA^+VII <C(||A^X/|| + ||A''/||) , (3.4) 







5(72 X 




T 


1)(2/-V, 



holds for every / in the Schwartz space. Looking at ( I3.4> with 7 0, we see that 

5 II II 2 

K has compact resolvent. Since e~2lli'll is an eigenfunction with eigenvalue for 
K*, it follows that K has also an eigenfunction with eigenvalue 0, let us call it g. It 
follows from (13 .4> and a simple approximation argument that || A''g|| < 00 for every 
7 > 0, and therefore g belongs to the Schwartz space. Furthermore, an argument 
given for example in I EPRB99. Prop 3.6] shows that g must be positive. Since one has 
furthermore 

the function p given by ( 13. 2> is the density of the invariant measure of ( 12. U . This 
concludes the proof of Theorem l3.ll □ 
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Before we give bounds on (12. 4> . we show the following little lemma. 

Lemma 3.2 Let S e (0, 2a^^) and let K be as in ( 13. it . Then, the kernel of K is 
one-dimensional. 

Proof. Let g G kerii'. Then, by the same arguments as above, e"^"^"^^ is the density 
of an invariant signed measure for Vt- The ergodicity of Vt immediately implies g oc g. 

□ 

Now we are ready to prove estimates on the solution of the Poisson equation (13. 1> . 

Theorem 3.3 Let h e C°°(T" x R") with D^yh e L^(T' x W;e-^\\y\\'' dxdy) for 
every multiindex a and every e > 0. Assume further that J h(x , y) ii(dx dy) = 0, 
where is the invariant measure for \2.1\ . Then, there exists a function f such that 
\3.1\ holds. Moreover, for every 5 > Q, the function f satisfies 

f(x,y)^e^\\y\\"f(x,y), feS. (3.5) 

Furthermore, for every 6 £ (0, 2ct^^), / is unique (up to an additive constant) in 
i2(T« X W^,'e-^\\y\\^dxdy). 

Proof. By hypoellipticity, if there exists a distribution / such that (12. 4> holds, then / 
is actually a C°° function. 

We start with the proof of existence. Fix 5 E (0, 2a^^), consider the operator K* 
which is the adjoint of the operator K defined in ( 13. 3> . and define the function 

— A II II ^ 

u(x,y) — h(x,y) e ^HJ'ii . 

It is clear that if there exists / such that K*f — u, then / = e^H^H^/ is a solution to 
(13. H . Consider the operator K*K. By the considerations in the proof of Theorem l3.1l 
K*K has compact resolvent. Furthermore, the kernel of K*K is equal to the kernel 
of K, which in turn by Lemma IJ!2l is equal to the span of g. Define Ti, = (g)^ and 
define M to be the restriction of K*K to H. Since K*K has compact resolvent, it has 
a spectral gap and so M is invertible. Furthermore, since Cy — T~^^^vix) — r^^y, 
one checks easily that f E H, therefore / — KM~^u solves K* f = u and thus leads 
to a solution to ( 13. 1> . 

Since K* satisfies a similar bound to ( 13 .4> and since || A'^m|| < oo for every 7 > 0, 
the bound (13. 5> follows as in Theorem l3.ll The uniqueness of 11 in the class of functions 
under consideration follows immediately from Lemma l?!^ □ 

Remark 3.4 Note that the solution / of ( 13. U is probably not unique if we allow for 
functions that grow faster than e"^ H^H . 

Remark 3.5 The identity yCp — 0, where C is the formal adjoint of C, immediately 
yields that j y fj,{dx, dy) = ^/t J v(x) ij,{dx, dy). In particular, the assumption that 
the drift is centered implies that y is also centered. Moreover, y clearly satisifies the 
smoothness and fast decay assumptions of Theorem l3.3l Hence, the theorem applies 
to each component of equation i2.4l and we can conclude that there exists a unique 
smooth vector valued function $ which solves the cell problem and whose components 
satisfiy estimate i3.5\ . 
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3.2 Estimates on the particle velocity 

One has the following bound 

Lemma 3.6 There exists a constant 7 > such that 

Eexp(i||a-iy(i)||2) <exp(i||a-i2/(0)f +7i 
Eexp(^^ \\a-^y{s)fds) < e^^(h^\a-^y{Q)f + lt 
holds for any initial condition y(0) and every t > 0. 

Proof. Itos formula yields immediately the existence of a constant 7 such that 



^ f \\a-^yis)f ds + ^ f {a-^y{s),dP{s)) 
2t" Jo Jo 



The first bound follows by exponentiating both sides and taking expectations. The 
second bound follows in a similar way after dividing both sides by 2. □ 

This yields the following: 

Theorem 3.7 Let ijj : T" x M" ^ M be such that 

1 



sup 



'0(a;,y)exp(--||CT ^y|p 



4 

Then, there exist constants C.S > such that 



< 00 



E{^j(x(t),ym - / ^ix,y)tiidx,dy) < C exp{\\a-^y{0)\\^ - St) . (3.6) 



T"xl 



Proof. From the smoothing properties of the transition semigroup associated to ( 12. 2K 
combined with its controllability and the fact that is a Lyapunov function, one 
gets the existence of constants C and S' such that 



||Pt(a;,y;-)-MllTV<C(l + ||y||V 



(See e.g. IIMT93I for further details.). Here — mIItv denotes the total variation dis- 
tance between the measures /z and v. Cauchy-Schwarz furthermore yields the generic 
inequality 

J fdpi- j fdv\ < ^\\^i-i,\\jy J fHdfi + diy) (3.7) 

The bound (13. 6t immediately follows by combining Lemma lT6l with ( I3.7l l. □ 
We also have a much stronger bound on the supremum in time of the solution: 

Lemma 3.8 For every k > and every T > 0, there exist constants S,C > such 
that 

E sup exp((5|l?/(s)||2) < C'e-'-e'^ll«(°>ll' , 
holds for every e € [0, 1]. 
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Proof. Let y be the Ornstein-Uhlenbeck process defined by 

1 /■* - 



IT _ ^ 

Then (see e.g. IAdl90l ). there exists constants c\ and C2 such that 



P( sup |iy(t)|| > a) < cie- 
^te[5-,s+T] ' 

for every s > 0. This immediately yields 



P sup ||y(t)|| > A < ci£-^e 

^tG[0,T£-2] ^ 

which in turn implies that there exist constants C3 and C4 such that 

e( sup exp(c3||?;(<)|p)) < C4£"^ . 

^te[0,Te-2] / 

The claim follows immediately by choosing 5 — {c^k)/2 and by noticing that there 
exists a constant C4 such that ||y(s)|| < |ly(s)|| + ||y(0)|| + C4 for all s > almost 
surely. □ 

4 Proof of Theorem O 

In this section we prove Theorem l2.1l 

Proof. By Theoreml33lwe have $(y, z) e C°°(T" x M", M"), so we can apply the Ito 
formula to the function , ) to obtain: 

£ Jo e V" Jo 

£ Vt Jq £ Vt Jo 



where we defined /3^(i) = ^i) and we used ( 12. 4> to get the second line. We also 

interpret Vy$ as a linear map from M" into W\ Thus we have: 

11/"* 

£ Jo 

= ea:-e($(yf,xf)-$(y,x)) + ^ / V y<^{ylxl) dp'{s) 

VT 7o 

ea; + + M^(i) . (4.1) 



It follows from (13. 5t and Lemma l378l that. for every p > there exists a constant C 
such that 

E sup |/i(i)r < Ce'^ . 

tG[0,T] 



It is therefore sufficient to show that \2.1\ and (12. 8> hold with /^e replaced by the law of 
the martingale term . We first show that ( 12. 7> holds. This is equivalent to showing 
that, for every £ e B* one can construct a random variable such that 

E|B£ - £(M")|P < C(r"f , (4.2) 
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holds uniformly over < 1, and such the law of Bg is given by We therefore 
fix ^ S B* with ||£|j < 1, which we interpret as a R"-valued measure with total mass 
(i.e. the sum of the masses of each of its components) smaller than 1. We also use the 
notation = £([t, T]). 

Integrating by parts, we can write 











We now define on the interval [0, T] the M-valued martingale M| by 

VT Jo 

According to the Dambis-Dubins-Schwartz theorem IRY99I Thm 1.6] there exists a 
Brownian motion B such that Mfit) can be written as 



i{s),{Vy^®Vy^){vl,xl)i(s))dS 



MI{t) = B{{MI,MI)t)^B 



On the other hand, the measure is a centered Gaussian measure with variance 
{i(s), K?l{s)) ds, so we can choose Bi to be given by 



Bf = B« 



Bl 



l^B(^f {e(s),IC^£(s))ds^ 



We wiU actually show a stronger bound than ( 14. 2> . namely we will show that 

JP := E sup \Bl - MI(t)\P < Ce"P . 
te[0,T] 



(4.3) 



We use the Holder continuity of the Brownian motion B, together with the Cauchy- 
Schwarz inequality to derive the estimate 

JP<E(miPjB) sup f\e(s)J—{Vy'P(^Vy^)iylxl)-JC^)£(s))ds"'') 

^ 0<t<T Jo \ ^ T J I J 



< (EHol2P(S))5(E sup / {l(s),H{yl,x\)i{s))ds 



< c 



(e su 

^ 0<t 



sup 



* 2aps i 

{l{s),H{yl,x\)l{,s))ds 







(4.4) 



where we introduced the n x n-matrix valued function 



H{x, y) = — (Vy$ ® Vj,$)(y, x) - JC'' . 

T 

In deriving the above estimate, we have used the fact that for a < i, the a-Holder 
constant of a Brownian motion is uniformly bounded on every bounded interval IRY99I 
Thm 2.1]. 

Note now that since £(t) is of bounded variation, £(t) (X) £(t) is also of bounded 
variation, so there exists a n x n-matrix valued measure £ on [0, T] such that £(t) (g) 
£(t) = £([t, T]). Therefore, we can integrate by parts in ( I4.4> to obtain 



JP < c(e sup Tr / / 

^ 0<t<T Jo Jo 



ft ns 2ap\ - 

H{y^^,x'^)dr£(ds) 
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<C E sup / H{ylxl)ds 

^ 0<t<T Jo ' 



Consider now the Poisson equation 



(4.5) 



By the definition of K?, we have / Fl(x, y) fj,{dx, dy) — (for each component), and 
we furthermore have exp(— (5||y|p)-ff € S for every 6 > 0. Therefore, using the same 
reasoning as in the proof of Theorem l3.3l equation ( 14. 5> has a unique smooth solution 
satisfying 

]\\vf Ti-r ni\ -T <= (4.6) 



F{x,y) ^ e--\\y\\ Hx,y) , TeS 
for every 5 > Q. We can apply Ito formula to deduce as before that 

/ H(ylxl)ds^~eHF{ylxl)-F{y,x)) + ^ [ VyF(yl,xl)a 
Jo Jo 



d(3(s) 



Therefore: 



|jp|2<^4apE \\F(yt,xt) 

te[0,T] 







te[0,T] 

Combining Lemma l3^ with ( 14. 6> . the first term can be bounded by 



yyF{yl,xl)d(3{s) 



2ap 



e^'^PE sup \\F{yl^t> 
te[0,T] 



In order to control the second term, we use the Burkholder-D avis-Gundy inequality 
followed by Holder's inequality, assuming that P > 



E sup 

te[0,T] 



VyF{ylxl)dP{s) 



2ap 



< CE 

^Ji 



\\\/yFiyl,x'M'ds 



sup E\\VyF(yt,xl)\\^''P 
te[0,T] 



This is bounded independently of e by ( 14. 6> and Lemma 13761 and so Jp < Ce'^P, for 
p > i. When p < one can bound Jp using the higher order moments. This 
completes the proof of bound ( I4.2l i and thus of the first part of Theorem l2.1l 

The proof of the second part of Theorem l2.1l is obtained in a straightforward way 
as a particular case of ( 14. 3> if one makes the choice I — kSx- □ 



5 Conclusions 

The problem of homogenization for periodic hypoelliptic diffusions was studied in 
this paper It was proved that the rescaled particle position converges to a Brownian 
motion with a covariance matrix which can be computed in terms of the solution of the 
Poisson equation (12. 4> . Further, an upper bound on the convergence rate in a suitable 
norm was obtained. Our analysis is purely probabilistic and this enables us to obtain 
more detailed information than what one could obtain from studying the problem at the 
level of the Kolmogorov equation. 

A very interesting question is whether a homogenization theorem of the form 12. II 
holds for random drifts v(x, t) and, if yes, under what conditions on v(x, t). From a 
mathematical point of view, it would be interesting to know whether it is possible to 
achieve convergence in the p-Wasserstein distance for n > 1. We plan to come back to 
these issues in a future publication. 
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